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ABSTRACT 

We show examples which reveal influences of spatial topologies to dynamics, using a class of 
spatially c/osed inhomogeneous cosmological models. The models, called the locally U(l)x U(l) 
symmetric models (or the generalized Gowdy models), are characterized by the existence of two 
commuting spatial local Killing vectors. For systematic investigations we first present a clas- 
sification of possible spatial topologies in this class. We stress the significance of the locally 
homogeneous limits (i.e., the Bianchi types or the 'geometric structures') of the models. In 
particular, we show a method of reduction to the natural reduced manifold, and analyze the 
equivalences at the reduced level of the models as dynamical models. Based on these fundamen- 
tals, we examine the influence of spatial topologies on dynamics by obtaining translation and 
reflection operators which commute with the dynamical flow in the phase space. 



^JSPS Research Fellow. Electronic mail: tanimoto@yukawa.kyoto-u.ac.jp 



1 Introduction 



Inhomogeneous spacetime models generally inherit from the full relativistic dynamics the strong 
nonlinearity in straightest ways, so that a moderately simple inhomogeneous model can be 
a good testing ground to obtain an insight toward understandings of generic properties of the 
relativistic dynamics. One of the best-known such models is the Gowdy model, where, first of all, 
two commuting spatial Killing vectors are assumed. By this assumption the spatial manifold 
is simplified to a one dimensional reduced manifold. Second, it is assumed that the spatial 
manifold is compact without boundary, i.e., closed. This assumption is favorable in view of the 
fact that no ambiguities occur in the boundary conditions once a topology is specified. A closed 
space is physically natural, as well, due to, e.g., the finiteness of the gravitational action. 

A striking property resulting from these two assumptions is the fact that the possible spatial 
topologies are restricted. When Gowdy models first appeared in 1971 [Q] as global generalizations 
of cylindrical plane wave solutions, only two spatial topologies, S'^ x and S^, were considered. 
After a few years from the first paper, Gowdy became aware of the fact that any closed three 
dimensional Riemannian manifold which admits two commuting Killing vectors is homeomorphic 
to one of T^, S'^ x S^, S"^ or a lens space L{p,q). So, he added model to his consideration. 
(Since lens spaces are finitely covered by the sphere S^, he argued little about lens space models.) 
This class of solutions or cosmological models consists of the original version of the Gowdy 
models. 

Recently, it has been pointed out p,^] that a natural generalization is possible, where the 
two commuting Killing vectors are local, i.e., the Killing vectors are defined in a neighborhood 
of every point but are not necessarily globally defined. (More precise definition of local Killing 
vectors will be presented in the next section.) This generalization gives us a set of new favorable 
models to test the dynamical properties of relativity, especially, as we will see, in connection 
with the spatial topology. 

Part of the motivation of this work comes from a plan to systematically investigate how 
the spatial topology influence the dynamics of a spacetime. In fact, the relativistic dynamics 
of a spacetime often seems to be influenced by its spatial topology. For example, we know the 
recollapsing conjecture [^^, where it is claimed that the well-known recollapsing property of a 
positive curvature (topologically S*^) homogeneous and isotropic cosmological model continues 
to hold also when the spacetime is inhomogeneous (if only appropriate energy and pressure 
conditions are fulfilled) . Thus we may interpret the recollapsing property is a direct consequence 
derived from the spatial topology. 

As another example, note (e.g., [^) that the dynamics of a vacuum homogeneous universe 
can be well characterized by the "potential" (given by the spatial scalar curvature), which can 
be determined if we know the Bianchi type. The local dynamics of a compactified (locally) 
homogeneous model can also be understood in the same manner (see |9-11]). Moreover in most 
cases, the corresponding Bianchi type is unique if the spatial manifold is closed. Thus we can 
say that the spatial topology determines the dynamical behavior of the universe at least at the 
locally homogeneous limit. 

How about when the spacetime becomes inhomogeneous? Can one, as in the locally homo- 
geneous cases, characterize the dynamics of an inhomogeneous universe by the spatial topology? 
When the spacetimes become inhomogeneous distinctions in local properties will obviously dis- 
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appear, as the Einstein field equations can be written in the same form, irrespective of the spatial 
topology. However, the boundary conditions may be different in general. So, we can naturally 
expect that the distinctions in the dynamical properties due to the differences in the spatial 
topologies manifest themselves in their global, discrete, statistical, or average properties. (Note 
that the recollapsing property mentioned above itself is also regarded as a global property.) It 
seems fascinating to find such properties. Also, such an investigation might provide us with 
some clues to the topology of our universe. 

In this paper we show, after some preparations availing ourselves of the generalized Gowdy 
models, an example of such properties using those models. More precise contents and contribu- 
tions should be the following. 

We first reveal all the possible spatial topologies of the generalized Gowdy models and 
classify them. In particular, we categorize all the models into two kinds, then subdivide each 
of them into finite types each of which has a correspondence to a particular type of locally 
homogeneous manifolds. Although we show the best classification we have, which will require 
further considerations (in fact, the number of the possible topologies is infinite), the set of 
the (finite) types mentioned above will be found to present the most relevant set of classes 
in an appropriate sense. (To be precise, however, some of the classes can have a discrete 
parameter.) While by our systematic analysis we find no new topologies other than those in the 
literature , one of our contributions is that we first asserts that those models do exhaust all 
the possible ones. 

Then we will concentrate on the models belonging to one kind, called the first kind, where 
each spatial manifold is a T^-bundle over S^. For these models we argue how we can represent 
the metric, and discuss the question of natural reduction (to as the spatial part). As a result 
of the reduction, we find that some models give rise to the same reduced Einstein equations 
with the same boundary conditions for their metric functions. We say that such two models are 
dynamically equivalent to each other. The dynamical equivalence greatly diminishes the number 
of representative models (or topologies) we should consider. The basic idea for the reductions 
is the same as the one presented in Ref. while the dynamical equivalence is first introduced 
in this paper. All of these settings are an indispensable step for the subsequent study of the 
models. 

Finally, motivated by the reason explained above, we consider varieties of spatial translation 
and reflection symmetries, and then ask if these symmetries imposed on initial data sets are 
preserved in time or not for each reduced model. In fact, we find remarkable differences in the 
reflections. 

These are done in Secs.^ to Section |5| is devoted to a summary and comments. In 
particular, comments on AVTD behavior and the local U(l) models are made. Appendix A gives 
a summary of the standard description of the Gowdy models with generalizations. Appendix B 
gives a summary of a calculation for the compactification of Sol geometry. 

We adopt the abstract index notation |jT^, that is, we use small Latin indices a, 6, • • • not 
to denote components but to represent the type of a tensor explicitly. To denote components 
we use Greek indices /i, • • • or capital Latin indices A,B,---. Conventionally we write tilde to 
denote a metric on a universal cover like gab , while a metric on a quotient space is simply written 
as gab- In the component representation, however, these metrics are possibly represented in the 
same way like g^u. We assume that all the spatial three- manifolds are orientable in this paper. 
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2 Topologies and Geometries of Locally U(l) x U(l) symmetric 
spaces and spacetimes 

Let (M, hab) be a Riemannian manifold . Suppose for any point p in M there exit neighborhood 
U which admits Kilhng vector fields, but these Killing vectors are not necessarily defined on 
whole M. For example, consider the manifold with metric 



where a{x) = a{x + 27r) and P{x) = f3{x + 2tt) are real periodic functions. This Riemannian 
manifold, denoted as M\ hereafter, possesses three independent Killing vectors 
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ix,y,z) {x + 2TT,-y,-z), 
{x,y,z) {x,y + 27:,z), 
{x,y,z) {x,y,z + 27r), 



(3) 



where = ^ is a vector field defined for convenience, we obtain a closed manifold Mi home- 
omorphic to T'^/Z2. (Here, an exponential of a vector represents the diffeomorphism generated 
by the vector.) See Fig.^. Note that ,^2 (and .^3) on the bottom points the opposite direction 
to that on the top, which fact tells us that ,^2 and .^3 cannot be defined on the whole Mi. ^4 is 
also incompatible with the identifications by the translations by 52 and (73. Thus, Mi admits no 
Killing vectors, though we can define Killing vectors on a neighborhood of every point. 




Figure 1: T^/Z2, generated from R'^ endowed with the Riemannian metric (P, by three identification 
generators gi, (72, and 173. 52 and gs are translations along, respectively, y and z-axes. gi is the composition 
of a translation along x-axis and the rotation in the y-z planes by tt. 

If we consider a spacetime manifold (Mi x R, gat) whose spatial metric coincides with Eq.(0), 
this gives a simplest (but rather uninteresting) example of a generalized Gowdy model. 

Definition Suppose a Riemannian manifold {M,hab) possesses the following two properties: 
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(i) There exists an open cover {Oj} of M such that every {Oi,h/^l) admits m independent 

Kihing vectors • • • , where h\l is the natural metric on Oi inherited from hah- 

(ii) On each Oi n Oj{^ 0), and ^(j) are related through a linear transformation 

m 

= « = l~n, (4) 

/3=1 

where /a"'^^ is a nondegenerate constant matrix (for fixed i and j). 

By local Killing vectors on (M, hab), we mean the collection of the pairs {{Oi], {^i^* ■ ■ ■ }) of 
the cover and the set of Killing vectors on each patch. For simplicity, we also simply denote 
them as ^1 • • • if we do not have to say on which patch these Killing vectors are defined. 

Prom the property (ii), the Killing orbits are well defined globally, even though the Killing 
vectors are local. 

Definition A three-dimensional closed Riemannian manifold (M, hah) is called a locally U(l) x U(l) 
symmetric space (or an CU'^ -symmetric space in short) if (M, hah) admits two commuting local 

Killing vectors [Jj({Oi}, l^^'-*, ^2'"*}); i-^-; [^i*'')?2*''] = Oon every Oj. We also impose the condition 
that this manifold must have no other local Killing vectors than those mentioned. We call this 
condition the genericity condition. 

As we will see later, the Killing orbits for this kind of space are flat 2-tori if the local Killing 
vectors are nondegenerate. We can choose each patch O, in the above definition as a regular 
neighborhood of such a (i.e., a neighborhood such that the boundaries coincides with another 
orbits), since the two commuting Killing vectors can be defined on such a neighborhood. In this 
case, the group U(l) x U(l) acts on this patch as isometries, hence the word "locally U(l) x U(l) 
symmetric" for the whole manifold (M, hah)- However, this terminology is rather lengthy to use 
frequently, so we also use 'XZ//^-symmetric" or more simply in this paper. 

Definition A spacetime manifold (M x R, gah) is called a locally U(l) x U(l) symmetric space- 
time (or a generalized Gowdy spacetime) if the spatial manifold M is closed and the spacetime 
admits two commuting spatial local Killing vectors. The genericity condition is understood, as 
in the case of the >CW^-symmetric spaces. 

Our first task is to determine all the possible topologies for the locally U(l) x U(l) symmetric 

spaces. It is convenient to consider separately the case where the local Killing vectors degenerate 
on some points and the case where they do not on any points. We call the latter (nondegenerate) 
type of manifolds the first kind, and the former (degenerate) ones the second kind. 

2.1 The first kind 

Lemma 1 // (M, hab) is an orientable locally U(l) x U(l) symmetric space of the first kind, 
then all Killing orbits must be closed and homeomorphic to T^. 
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Proof. Let (O, rjab) be a Killing orbit, where r]ab is the induced metric on the orbit from 
hab- Since the Killing orbit possesses two commuting Killing vectors, rjab is a flat metric. Since 
M is orientable, O is homeomorphic to one of R^, R x 5*^, or T^. First, let us assume O ~ R^. 
Consider a geodesic / in (O, r]ab), and consider a point sequence {qi} such that all are on I and 
and qi^i are unit distant with respect to rjab- In M the sequence {gj} must have a converging 
subsequence {pj}, since M is compact. Let be the neighborhood of the limit point poo with 
radius e, i.e., |p — poo\ < e for all p G where | • | is the geodesic distance with respect to hab- 
Ve contains all pi {i > N^) for an integer N^, so that for any pi and pj, > N^), \pi — Pj\ < 2e 
as a result of the triangle inequality. On the other hand, the geodesic / is a "straight line" in 
the Euclid plane (R^,ryaf))5 so that the distant between any two points pi and pj (i / j) in O 
is larger than or equal to unity. Hence we can choose a sequence of "isometric disconnected 
neighborhoods" in M such that pi G Ui, Si n Sj = 0, where Si = UiCi O, and there exist 
an isometry cpij : {pi,Si,Ui) {pj,Sj,Uj) for all i,j. Since Si Ci Sj = 9 and \pi — Pj\ < 2e, 
we must have the consequence pj Si but pj G Ui for a sufficiently small e (i.e., for sufficiently 
large i and j). Moreover, this implies that for an arbitrarily small e' < e there exist two points 
in M, Pi and pj, such that \pi — Pj\ < 2e', pj 5^ but pj G ?7j, and there exist the isometry 
<Pij ■ Pi Pj- This in turn implies that there exist a third motion in {M,hab) off the orbit, 
which fact is against the genericity condition. (That is, R^-orbits can be formed only if (M, hab) 
is locally homogeneous.) In the case of O ~ R x S^, we can choose a geodesic / which is open as 
in the case of O ~ R'^. Repeating a similar argument we conclude that this case is also against 
the genericity. Thus, the only possibility is the case O ~ T^. □ 

Theorem 2 // {M,hab) is an orientable locally U(l) x U(l) symmetric space of the first kind, 
then M is homeomorphic to a T"^ - bundle over S^ , for which the local Killing vectors generate 
the T"^ -fibers- 
Proof. Consider a normal vector field that is everywhere nonvanishing and normal to the 
Killing orbit, and consider the integration curve c of this field passing through an arbitrarily 
given point p G M. Let O be the Killing orbit to which p belongs, so c intersects with O at p. 
If c did not intersect again with O, then it would imply that the total volume of M is infinite, 
which fact contradicts with the compactness of M. Hence c intersects with O again. This 
implies that if we regard each Killing orbit as a point, then M reduces to S^. Since all orbits 
are homeomorphic to as in Lemma |l], M is a T^-bundle over S . □ 

Topologically, any T^-bundle over S^ can be obtained by first considering the product x /, 
where / = [0, 1] is the unit interval, then identifying x {0} and x {1} by the action of 
a gluing map (p : T'^ ^ T'^. Since any two gluing maps which are homotopic to each other 
results in topologically the same manifold, we usually think of the gluing map (j) as an element 
of the mapping class group of T^, mcg(T^). Here, the mapping class group of a manifold M 
is the group of diffeomorphisms of M modulo diffeomorphisms which are homotopic to the 
identity, mcg(M) = Diff (M)/Diffo(M). The group mcg(T^) is well known as the modular 
group, GL(2, Z). However, since we are interested only in orientable manifolds, we consider the 
orientation-preserving mapping class group of the torus, mcg^(T^) ~ SL(2, Z). Letting 




(5) 
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the fundamental group of the corresponding space is represented by 



TTi = (91,52,53; [92,93] = 1, 515251 ^ = 5253, 5i535i ^ = 5253/ , (6) 

where 52 ™d 53 are generators of the fiber and gi is that of the base S^. We therefore have 
a natural map loi : SL(2, Z) — > Wi, where Wi is the set of orientable T^-bundles over S^. While 
cji is not injective, we have the following useful theorem, of which proof can be found in Ref. [ p^ 
(Theorem 5.5 therein); 

Theorem 3 ( [|13|) Let M be the total space of a T"^ - bundle over with gluing map <f), and 
let A G GL(2, Z) represent the automorphism of tti{T'^) induced by 4>. Then, M admits a Sol- 
structure if A is hyperbolic, an E"^ -structure if A is periodic, or a Nil-structure otherwise. In 
particular, i/|Tr^| > 2, then A is hyperbolic, so M admits a Sol- structure. 



This theorem tells us that every T^-bundle over admits a "geometric structure" in 
other words, admits a locally homogeneous metric. This structure is, as in the above theorem, 
one of the three types, E^, Nil, or Sol. Conversely, we can find all the T'^-bundles over from the 
closed quotients oi E^, Nil, and Sol, so we can refer to known classifications of the closed quotients 
of the three geometries to classify T^-bundles over S^. This procedure simultaneously determines 
the corresponding geometric structure for each representative. (To avoid confusions we remark 
that locally homogeneous metrics are used just for convenience here to discuss topologies, but 
we will see that appropriate inhomogeneous metrics are obtained by "relaxing" such a metric.) 

The case of E^: All the orientable closed manifolds modeled on E^ are classified into six 
manifolds H, T^, TVZ2, TVZ3, TVZ4, TVZe, and T^ /Z2 x Z2. All of these except the last 
one can be regarded as T^-bundle over S^. Associated with these five are the representative 
elements Ex £ SL(2,Z) (A = 1, 2, 3, 4, 6, respectively): 




^1 = n ; , ^2 = n \ \, E: 





(7) 

The case of Nil: All the orientable closed manifolds modeled on Nil can also be completely 
classified (See page 4878 of Ref. P). Among them, ones which can be regarded as T^-bundle 
over are given by the following two families A'^-i-i(n) G SL(2,Z): 





where n is a nonzero integer. The parameter n in Ni(n) can be chosen to be positive, since 
the relabeling (91,93) (53,51) reverses the sign of n. For A^_i(n), the opposite sign of n 
corresponds to a distinct topology. 

The case of Sol: As in Theorem ^, if \TtA\ > 2, we have a T^-bundle over modeled on 
Sol. All those for \TvA\ < 2 are non-orientable so we will not consider them. Any two ones 
with distinct values of Tr^ such that \TrA\ > 2 are not homeomorphic, so one parameter family 
0tSL(2,Z) @, 
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where n is an integer such that \n\ > 2, gives a one parameter family of distinct T^-bundles 
over modeled on Sol. However, this does not implies that the ones modeled on Sol can be 
classified only by TiA like Eq. (^) , since there exist topologically distinct manifolds with the same 
TV A [16|. No complete classification of the closed manifolds modeled on Sol seems to be known. 
Nevertheless, as we will see in the next section, the sequence (|9|) gives nice representatives as 
reduced relativistic models. 



2.2 The second kind 

Theorem 4 If (M, hab) is a locally U(l) x U(l) symmetric space of the second kind, then M is 
homeomorphic to one of S^, 5^ x , or a lens space L{p,q). 

Proof. Since M is closed, if removing an open neighborhood of the degenerate points from 
M, the resulting manifold M with boundaries is compact. Moreover, the boundaries can be 
chosen so as to coincide with (regular) Killing orbits. Let dM be the boundaries of M so 
chosen. Note that in the proof of Lemma |l|, only the compactness of M is assumed to show 
that the Killing orbits are T^. Hence all the Killing orbits for (M, hab) are also T^. (Here, 
hab is the restriction of hab to M.) In particular, every connected component of dM is T^. 
Consider a (nonvanishing) normal vector field with respect to the Killing orbits in (M, hab)- This 
flow of the normal vector field defines a unique one-to-one correspondence between boundary 
components, since the image of a boundary component by this flow must end up with another 
boundary component, due to the compactness of M. Since M and therefore M are assumed 
to be connected, M must have only two boundary components and is naturally homeomorphic 
to the product x /, where I = [0, 1] is the unit interval. Next, consider the neighborhood 
of a degenerate Killing orbit removed. An action of U(l) x U(l) can degenerate only to U(l), 
so a connected set of degenerate points forms a circle, and a neighborhood of such a circle 
with boundary forms a solid torus. The boundary, which is homeomorphic to T^, can again be 
chosen so as to coincide with a (regular) Killing orbit. Consider two such neighborhoods, Vi 
and V2- Then, the original manifold is recovered by attaching them to M along the boundaries: 
M ~ Vi U M U 1^2. Note, however, that Vi U M is another sohd torus with the U(l) x U(l) 
symmetry. Let us rewrite this manifold as Vi. The original manifold M is now obtained simply 
by identifying the boundaries of the two solid tori Vi and V2, M ~ U V2- Finally, it is a 
well-known fact (e.g. ||T^) that the sum of two solid tori is homeomorphic to one of 5"^, 5"^ x 5^, 
or a lens space L{p,q). □ 

Cosmological models based on £Z//^-symmetric spatial manifolds of the second kind can be 
thought of as a global generalization of cylindrically symmetric (plane wave) models. In fact, the 
symmetry axis of a cylindrically symmetric space is a degenerate orbit for the U(l) x U(l) action 
of the cylindrical symmetry. If, as usual, taking this axis as z-coordinate axis, and introducing 
identifications z ^ z + zq for a constant zq, we have the symmetry axis that is a circle. A regular 
neighborhood of this circle is a solid torus with the U(l) x U(l) action. An £W^-manifold M 
can be obtained by identifying (the boundaries of) such two U(l) x U(l) symmetric solid tori. 

Now, it should be remarked that the local Killing vectors for an £Z//^-manifold of the second 
kind are always globally defined. This is because this kind of space can also be regarded as 
the product x / (with a metric that degenerates on the boundaries). The original Gowdy 
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models in fact include all the models based on S^, S'^ x S^, and the lens spaces L{p,q). We 
therefore do not discuss much about these models in subsequent sections. 

We end this subsection with a summary of the classification (e.g. [|l7|) of the lens spaces, 
for completeness. Consider a gluing map </> : ^ to identify the boundaries of two 
U(l) X U(l) symmetric solid tori, dVi and dV2. Since if two gluing maps are homotopic the 
resulting manifolds are homeomorphic, it is again natural to think of an identification map 

as an element of the mapping class group, G mcg(T^) ~ GL(2,Z). Let ^ q ^ ^^'^ ( ^ ) 
be, respectively, meridian and longitudinal loops of dVi or dV2- We define the action of j4 = 

(V T \ 
1 G GL(2,Z) by the left action to these vectors. Then, the closed manifold M made 

by gluing the two solid tori with respect to A is the lens space L(p,q). (The topology of M is 
regardless to r and s, i.e., it is determined only by the mapping of the meridian loop of dVi.) 
The integers p and q are coprime and we can set p > (since L{p,q) ~ L{—p,q)). and 
S*^ X correspond, respectively, to L{l,q) and -L(0, 1). Lens spaces are completely classified by 
the fact that L(p, q) and L{p' , q') are homeomorphic if and only \i p = p' , and q = zizq' mod p or 
qq' = ±1 mod p. 



3 Relaxation, Reduction and Dynamical Equivalence 

In this section we present an detailed account of the reduction procedure for the £Z//^-symmetric 
models of the first kind, with emphasis on the significance of the corresponding geometric struc- 
tures. We will also obtain useful results, the "dynamical equivalences". 

First, we remark that since the spacetime manifold is the product M x R, where the spatial 
manifold M is a (T^-)bundle over with fibers generated by the local Killing vectors, the 
natural reduced manifold obtained by contracting the Killing orbits should be x R. However, 
as we see bellow, when the bundle is not (covered by) a direct product, this reduction is not 
trivial. 

We represent an £Z//^-symmetric spacetime (of the first kind) (M x R, g^;,) with the covering 
map 

TT : (M X R, gab) ^ (M X R, gab) = (M X R, gab)/r, (10) 

where (M x R, g^b) is the universal covering manifold of (M x R, g^;,), and F is a covering 
group, which acts on M discretely. Since the spatial manifold M is a T^-bundle over S^, M is 
homeomorphic to R^. We can therefore use globally defined coordinates, say {x^ 
{t,x,y,z), to represent gab- A natural form of the metric admitting two commuting Killing 
vectors represented with these coordinates is g^uit,x)dx^dx'^ , where Greek indices iJ,,^, - ■ ■ run 
over to 3. Moreover, fixing the freedom of diffeomorphisms in ^ij^i,{t^ x^dx^ dx^ cis much as 
possible, we can represent the metric gab &s O] 

ds^ = e^/'^{-dt^ + dx^) + 2NAdx^dt + ReAB^x^^x^ , (11) 

where 7, R, Na and cab (Capital indices A,B,--- run 2 to 3) are functions of t and x. We set 
deters = 1 so that R describes the area of the Killing orbit. This metric is called the generic 
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metric of the canonical form. We will also consider the restricted metric with Na = (the 
"two-surface orthogonality" [^) 

= e^/2(_dt2 + ^^2^ + RcABdx'^dx^, (12) 

which is called the two-surface orthogonal metric of the canonical form. For a standard pre- 
scription of Einstein's equations for the last metric, see Appendix (Note that the metric 
functions are represented with "bars" in this Appendix. We will take this notation from the 
next subsection to avoid conflicts with "noncanonical" metrics which will appear.) 

Both metrics admit two dimensional isometry group consisting of translations along y and 
z axes. Note that the y-z planes descend to the T^-fibers after appropriate identifications in 
each y-z plane. These identifications are generated by two independent vectors, which we can 
choose without loss of generality as the unit coordinate generators, d/dy and d/dz. Then we 

can naturally think of the actions of these generators as the representation of the generators, 

A. d_ 

§2 and 33, of TTi into the isometry group; §2 — e^y , 53 ~ es^. The remaining generator gi 
must induce a translation along x-axis plus the modular transformation induced by an element 
A = n^B £ SL(2,Z). Thus, gi : {x,y,z) {x -\- 27r,n^BX^) on the t =constant space. 
Together with (73 : (x, y, z) {x,y + 2it, z) and g^ : (x, y, z) (x, y, z + 27r), all the relations in 
the fundamental group (^ are fulfilled. The boundary conditions for the metric (11) or ( [l^ ) is 



determined by the requirement that the action of gi be an isometry of the metric, which can be 
easily found; 

X{t, x) = X{t, X + 27r), R{t, x) = R{t, x + 2it), 

NA{t,x) =n^ANB{t,x + 2TT), eAB{t,x) = rf au^ sec D{t,x + 2t:). (13) 

Note that the components cab that describe each T^-fiber are not periodic functions in general, 
so that the reduced manifold (spanned by t and x) cannot be naturally regarded as x 
R, if we represent the metric as Eq.(pT|) or (|l^). The reason why cab do not automatically 
become periodic is that what two independent components in cab themselves describe are 
Teichmiiller parameters for the T-^-fiber, rather than moduli parameters. (For the difference 
between Teichmiiller and moduli parameters, see, e.g., Ref. |l^.) 

An idea to remedy the situation is to choose the metric functions so that they become 
constant when the spatial metric is at a locally homogeneous limit. That is, conversely, we 
"relax" a locally homogeneous metric in an appropriate way to obtain a suitable metric. Note 
that a constant function is naturally a function on . As far as considering smooth deformations 
of the metric, metric functions chosen in such a way must continue to be (smooth) functions on 

even when the metric becomes inhomogeneous (with the CU^ symmetry). 

The scheme One nice way to realize this idea is to expand the spatial metric in terms of 



the (left) invariant one-forms a* (z = 1 ~ 3) of a Bianchi group G. As we have seen in Sec. 2.1, 
any T^-bundle over admits one of , Nil, or Sol-structure. These geometric structures 
correspond, respectively, to Bianchi I(VIIo), II, and VIq- {E^ has multiple correspondences.) 
The appropriate Bianchi group G is determined from this correspondence. Let (i = 1 ~ 3) 
be independent generators of G, that is, are Killing vectors of a G-invariant metric. Since 
(T* are by definition invariant under the action of G, the homogeneous metric is written as 
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hija^a^ with (spatially) constant components hij. The Killing vectors of this homogeneous 
metric include a commuting pair. Suppose [^2, '^3] = 0. Let x be a vector field such that x 
is independent of ^2 and ^3 at every point and is invariant under the action generated by ^2 
and ^3, i.e., C^2X = ^^aX = 0- Here, £ represents the Lie derivative. Since [^2,^3] = 0, we 
can choose coordinates {y and z) generated by ^2 and ^3. Moreover, since ^^.^x = ^^sX = 
implies the commutativity [C2,x] = [Cs^x] = 0) we can choose a third coordinate (x) as that 
of generated by x- Consider an arbitrary function f{x) that depends only on x chosen in this 
way. Then, the level set of f{x) is invariant under the diffeomorphisms generated by ^2 and 
1^3, and it coincides with the set of the orbits generated by .^2 and ^3. Hence we can write the 
inhomogeneous metric as hij{x)a^a^ , which is invariant under the diffeomorphisms generated by 
^2 and ^3, and inhomogeneous in the desired manner. (Recall that all are invariant under ^2 
and ^3.) The metric functions hij{x) should be periodic in x because of the reason explained 
above. Hence, a spacetime metric with hij{t, x)a^a^ as the spatial part naturally gives rise to a 
reduction onto x R. We call this scheme the relaxation method. If ^1 = x as in the Bianchi I 
case, then this scheme descends to the usual coordinate representation like Eq.([Tl|) or ([T^), but 
otherwise it is nontrivial. □ 

Recall that a Bianchi group G is a three-dimensional simply transitive group acting on 
a three-dimensional simply connected manifold M, and therefore we can identify G with M 
(e.g. ||l2[). We adhere to this viewpoint in this paper and use the components of G to represent 
the coordinates of M, as well. 

One comment should be made here. As pointed out in Ref. j^, all Bianchi homogeneous 
spaces except for VHI and IX possess a commuting pair of Killing vectors. Even for VIII and 
IX, if considering higher symmetry there can exist commuting Killing vectors. Moreover, these 
homogeneous spaces can be compactified to closed spaces, except for IV and Via types. Therefore 
all the locally homogeneous manifolds of the Bianchi types except IV and VI^ seems to be locally 
homogeneous limits of £Z//^-manifolds. However, as we mentioned above, only Bianchi I, VIIq, 
II, and VIq types can be such a limit for the first kind manifolds. The second kind manifolds 
correspond to Bianchi IX (in case of M ~ S*^ and L{p, q)) and the Nariai-Kantowski-Sachs type 
(in case of Af ~ S"^ x S*^). Thus, Bianchi III, V, Vila, VIII types are missing in the list of possible 
manifolds as the locally homogeneous limits of our CU"^ -manifolds. The reason is because each 
orbit of the two commuting local Killing vectors for a locally homogeneous manifold of these 
types always do not close. For example, in the case of Bianchi V, it forms R^. As we have seen 
in the proof of Lemma [l|, such a case is possible only when the manifold is locally homogeneous. 
That is, we can say that, as a corollary to Theorems |2| and 

Corollary 5 Every closed locally homogeneous manifold of Bianchi III, V, Vila and VIII types 
cannot be relaxed to be inhomogeneous with CU^ symmetry. 

In this sense, these manifolds are "stiff." It is worth noting that a common feature of them is 
that they all contain a hyperbolic structure, or . 

In the following subsections, we perform the reduction procedure for each type of Nil, Sol, and 
E^. The spacetime manifold M x R reduces to S"^ x R for any spatial manifold M. One important 
feature that is found as a result of the reductions is that, even for models with topologically 
distinct spatial manifolds, distinctions at the reduced level can completely degenerate. In other 
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words, we obtain for some models the same set of PDEs (as the reduced Einstein equations) 
with the same boundary conditions. For such models the universal covering spacetime metrics 
can be represented in exactly the same form with only the difference contained in the covering 
group r. We call that such models are dynamically equivalent to each other. 



3.1 Case of Nil 



The Bianchi II group Gn is the three-dimensional simply transitive group with the multiplication 
rule 

a a + X 



y 

z 



h + y 
y c + z + ay 



This group is generated by 



d d , 
ox oz 



d_ 

dy 



, ^3 



d_ 

dz 



The one-forms which are invariant under the action of Gu are 



a 



dx, a = dy, a = dz — xdy. 



(14) 



(15) 



(16) 



Since ^2 and ^3 commute and they generate coordinates y and z, a desired representation of the 
spatial metric is given by hij{x)cr^a^ . 

Before proceeding further, however, we have to clarify a subtle point about the fact that we 
can choose another commuting pairs of group generators. For example, ^1 and ^3 also commute 
to each other, and the orbits they generate are distinct from those generated by ^2 and ^3. In 
general, since the commutation rule for linear combinations of the group generators is given by 



(17) 



any two generators a = a*^i and (3 = P^^j commute if and only if a^P"^ — a^fS^ = 0. The 
last condition implies that ^3 must be tangent to the surface spanned by a and /? if a and /? 
are independent. Such a surface, conversely, is spanned by r]2 = sin^^i + cos and 773 = ^3, 
where & is a real parameter. Hence, we have one-parameter (6) family of distinct sets of orbits 
generated by a commuting pair. This then implies that we have one-parameter degree of freedom 
of relaxing a locally homogeneous manifold of Bianchi II type. 

However, this fact is actually insignificant if viewing the metric as a universal cover metric. 
That is, the freedom of 9 can be absorbed in the freedom of choosing the covering group P. 
Therefore we do not have to consider the freedom of 9, which can be fixed 9 = without loss 
of generality. (The model based on 9 = 7r/2 was presented in Ref. as 'Type 2', which is 
redundant for this reason.) 

Now, we have established the fact that the relaxed (universal covering) metric can be rep- 
resented by hij{x)a^a^ with the basis ( p!6[). The form of spacetime metric corresponding to the 
two-surface orthogonal class of metrics (|85| ) is given by 



2\2i 



(18) 
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where the metric functions 7, R, P and Q are functions of t and x, and periodic with respect 
to X. (We have changed the choice of P and Q from that in Ref. [P, since the present choice 
gives us the most natural and simplest correspondence to P and Q. See bellow.) These metric 
functions are related to those for the canonical representation (|85| ) through 



^ = ^^ R = R^ P = p^ Q = Q 



(19) 



The reduced Einstein equations for the unbarred variables 7, P, Q, and R are therefore easily 
obtained by substituting these equations into those for the canonical (i.e. barred) variables 
presented in Appendix A. 

Our next task is to check that we can indeed compactify the universal cover. We choose the 
period for the metric functions along the x axis as 2ti: 

f{t,x) = f{t,x + 2TT), for / = 7, i?, P and Q. (20) 

We denote the isometry group for the spacetime metric (p^) as -f^n ^ , and its subgroup which is 
a subgroup of Gu as Hu C Gu. We have 

/ 27rn 
b 



(2) 

The full isometry group Hj^ is generated by Hu and the Z2-isometry 

k : {x,y,z) ^ {x,-y,-z). (22) 
Recall that the fundamental groups for Ni{n) and A^_i(n) are given by Eqs.(P) with (0). Putting 



V 



n e Z,6,c G R 



(21) 



9i 



(23) 



we have to solve the relations in these fundamental groups for the parameters gi^ and n^. The 
solutions are easily found for Ni{n). For A^_i(n) we have to think of gi as a composite of k and 
an element of Hu- The solutions are given by 
For Ni{n): 



Tn = {91,92,93} 





(24) 



where p,qGZ, gi^ £ R, and pgs'^ - qgi'^ / 0; 
For N-i{n): 




(25) 
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where p G Z, gi^ G R, and pgs^ ^ 0. The point here is that there exist solutions, ensuring that 
we can compactify the universal cover for any closed spatial manifold N±i{n). For completeness 
we remark that the parameters appearing in r„ can be fixed arbitrarily, since redefinitions of 
metric functions can make the parameters equal to arbitrary values. For example, we can take 
p = g^'^ = 1, and q = g^^ = g^^ = g^^ = Q for iVi(n), and p = g^'^ = 1, and gi'^ = gi^ = g^^ = 
for iV_i(n). 

Note that we did not have to impose any other condition than Eq.(20) for all spatial topologies 
N±i{n). This results in obtaining the same reduced Einstein equations with the same boundary 



conditions pQ). Thus: 



Proposition 6 All (two-surface orthogonal) CU^ -symmetric models of Nil type are dynamically 
equivalent. 

Here, the two-surface orthogonality is necessary for N^i{n) models to ensure that the map k is 
an isometry. That is, generic A^_i(n) models with shift functions are not allowed. On the other 
hand, all generic Ni{n) models are dynamically equivalent. 

3.2 Case of Sol 

The invariant one-forms are given by 

cji = dx, = ^(e^^dy + e-<^^dz), = ^(-e^^dy + 6"^^^), (26) 
V 2 v2 

where q > is a positive parameter introduced for convenience. 

Using these 1-forms we can write the two-surface orthogonal metric as 

dc,2 = e^^H-dt' + {a')^) + i?[e^(a2' + Qa^'f + e-^'ia^'f], (27) 

where 7, P, Q, and R, are functions of t and x, and they are assumed to be periodic in x. We 
have defined 

a^' = -^{a^- a^) = e^My, a^' = ^{a^ + a^) = e'^^dz. (28) 
v 2 \/2 

(The choice of metric functions P and Q here is different from that in Ref . . We have done so 
for apparent simplicity, but we should keep in mind that at the homogeneous limit the metric 
is not diagonalizable with respect to {a^,a'^',a'^'), in contrast to with respect to {a^,a'^,a^), 
for vacuum spacetimes.) The metric functions defined with Eq.(p7|) are related to those in the 
canonical representation (|8^) through 



7 = 7, R = R, P = P + 2qx, Q = e-^'^^Q. (29) 

The reduced Einstein equations for 7, P, Q, and R can be obtained by substituting these 
equations into those presented in Appendix 

For compactifications, see Appendix ^ In particular, the isometry group of the metric (|27| ) 



(2) 

is given by the -ffyio Pi'ssented in the Appendix. Therefore we can compactify the universal cover 
possessing this metric for all A. If we choose the period along x-axis for the metric functions as 
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2tt, we must put C3 = 2tt. Accordingly, we must choose the parameter q so that e^'^^ equals to 
the greater absolute eigenvalue of the matrix A, since we have assumed q > 0. Note that the 



characteristic polynomial (107) depends only on TrA. Let n = Tr^. Then we obtain 



1 |n| + \/n2 -4 
q=^log . (30) 

The reduced Einstein equations for the (unbarred) metric functions depends (only) on the pa- 
rameter q, and the boundary conditions are the common simple periodic boundary condition. 
Moreover, q depends only on |n|. We thus conclude: 

Proposition 7 An CU^ -symmetric model of Sol type is specified with an element A G mcg_|_(T^) 2 
SL(2,Z) such that \TicA\ > 2. Let Ai and A2 he such matrices. Then, the corresponding models 
are dynamically equivalent if \ TiAi \ = \ Tr^2 1 • 

From this fact, we find it suffices to consider the models with the one-parameter family S{n) 
(n > 2), defined in Eq.(|9|), as representatives. 

3.3 Case of 

The corresponding Bianchi types are Bianchi I and VIIq. More precisely, geometry E'^ is the 
Euclid space Y{? with the standard metric dx^ + dy^ + dz^ . The isometry group of is therefore 
formed by the translations R'^ and rotations 0(3), and is isomorphic to the Poincare group, 
IsoTuE^ ~ 10(3). 10(3) contains two simply transitive subgroups, the Bianchi I {Gi ~ R^) and 
VIIq (GviIo) groups, so that the correspondence above is arrived at. 

We can make the inhomogeneous metrics from both Gi and GviIq- However, the effects 
are not equivalent. First, consider the metric made from Gj. Since the invariant one- forms 
coincide with the usual coordinate basis dx, dy, and dz, the spacetime metric and the boundary 



conditions are the same as the canonical ones (O) and (13). We find in particular that the 



metric functions for Ei and E2 are simply periodic. For £^3, £'4, and Eq, however, the metric 



functions must obey the last condition in Eqs.(13), which does not imply the simple periodic 



conditions, (though they have period Gvr, 47r, and Gvr for £3, £4, and Eq, respectively). As 
long as we adhere to the metric obtained from Gi, this feature is inevitable, since the closed 
manifolds corresponding to £3, i?4, and Eq cannot be realized in the Bianchi I group Gi, or 
Gj . Here, Gj is generated by Gj and the Z2-map h : {x,y,z) — > {x, —y, —z), which is an 



isometry of the metric (0). Precisely, the fundamental group 7ri(£'i) corresponding to Ei can 
be embedded in Gj, and vri(£'2) can be embedded in the Gj . However, other fundamental 

groups can be embedded neither in Gj nor Since the isometry group of the metric ([l^ ) 

(2) 

becomes a subgroup of Gj , provided that the metric functions are all periodic with period 27r, 
this implies that the embeddings of vri for £3, £4, and Eq into the subgroup are impossible. 
Thus, the appropriate boundary conditions cannot be the simple periodic conditions for them. 

On the other hand, all the fundamental groups -ki{Ex) (A = 1,2,3,4,6) can be embedded 
into Gviio [0) which fact implies that the boundary conditions can become simple periodic ones 
for the metric obtained from GviIq for ah Ex- We present an explicit prescription bellow. 
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The invariant one-forms for Bianchi VIIq are given by 

= dx, 0"^ = cos qx dy + sin qx dz, = — sin qx dy + cos qx dz, (31) 

where q > is a positive parameter. Using these 1-forms we can write the two-surface orthogonal 
metric as 

= e7/2(_dt2 + (^1)2) ^ ^[gP(^2 ^ g^3)2 ^ e-^^^a^f], (32) 

where 7, P, Q, and i?, are functions of t and x, and they are assumed to be periodic in x with 
period 2tt. The parameter q wih be chosen so that we can compactify the universal cover. (The 
role of q is similar to the q appearing in the Sol model.) The correspondence to the barred 
metric functions in the canonical metric (^) is given by 

^Pi^) = <I)(P(x), (5(2;), cos qx, sin qx), Q(x) = ^'(P(x), Q(x), cos qx, sinqx), 

gPi^) — <I)(p(x), Q(x), cos qx, — sinqx), Q{x) = ^'(P(x), (5(x), cos qx, — sinqx). (33) 

where the two functions $ and ^ are defined by, for arbitrary P, Q, c, and s, 

^>(P, Q, c, s) = e^(c - sQf + e-^s^ 

vi/^/^,C^,c,sJ $(P,Q,c,s) ■ ^ 

Note that the vacuum Einstein equations in terms of the unbarred variables contain the param- 
eter q. 

As usual, we represent the elements of the Bianchi group in the column vector form. The 
multiplication rule for GviIq is given by 

where Pqa is the rotation matrix by angle qa; 

P,.= f^°^^" -sinqa\ 
\ smqa cosqa / 

The subscripts q appearing in the column vectors are to remind that the multiplication rule is 
defined with respect to q. The one-forms (^) are invariant under this action. 

We denote the isometry group of the metric (|3^) which is a subgroup of GviIq as -ffyiio • It 
is given by 

27rn 







\ X 








y 




V 









Hyug = { \ b e Gviio 



c 

r(2> 



n E Z, b,cG'R} . (37) 



q 



The (full) isometry group -f^viio metric (^) is generated by ^^viio the Z2-isometry 

h : {x,y,z) {x,-y,-z). 
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We can perform the embedding of every 7ri{Ex) both with and without the Z2-isometry h, 
as we can read from Sec.V A of Ref. Q. (But, be careful with the differences in the choice 
of representations of the relations in tti.) When we do not use h, i.e., embed the fundamental 
group ■Ki{E\) into -ffyiio) the appropriate value of q is simply given by 

q=^. (38) 

Hence all the models seem to be dynamically inequivalent. However, we can choose q = 1 for 
A = 2, and q = 1/3 for A = 6, if we use h. For example, in the case of Eq one can confirm that 
the solution of the embedding into i^viio is given by 

Te = {91,92,93} 



















[ 










[si') 


: 

6 


6 







±7r/3 



92 
92' 



(39) 



with q = 1/6. Here, gi^ are real parameters. Also, we can embed the fundamental group into 



2tt 
2 



ho \ gi 

51' 



92 

92' 



R 



'2tv/3 



92 
92' 



(40) 



with q 



1/3. Since the fundamental group of E^ can be embedded into -fTviio with q ■ 

I 



3 as 







f 












92' 


91' ) 




{92' 



V 



-R 



±2n/3 



92 
92' 




(41) 



the Eq and E^ models are dynamically equivalent. (The significance of the parameters gi^ is the 
same as in the case of Nil. In particular, we can take gi^ = gi"^ = 0, (72^ = 1, and (72^ = 0, in all 
cases above.) 

As a result, we have only three dynamically inequivalent classes; 

Proposition 8 The Ei, E^, and E4 models comprise a set of representatives for the dynamical 
equivalence in the two-surface orthogonal CU^ -symmetric models of E^ type. 

For the generic models of E^ type with the shift functions, all the E\ models are dynamically 
inequivalent, since the map h is not an isometry. 

To end this subsection, we remark again that the £i^^-symmetric models of E^ type best 
match with Bianchi VIIq, not Bianchi I. However, if limited to the (or T^/Z2) model, which 
corresponds to Ei {E2), what it naturally corresponds is Bianchi I. To make later discussions 
easier, it may therefore be convenient to put forward the following convention. 

Convention When speaking of the models" , we understand them to have correspon- 
dence to Bianchi VIIq. In particular, the "unbarred variables" for an E"^ model are those in the 
metric (^). However, when speaking of the "T^ model", we understand it to have correspon- 
dence to Bianchi I. Therefore the unbarred variables coincide with barred variables appearing 
in the canonical metric (^). 
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4 Translation and Reflection symmetries 



The Gowdy equations (92) have several symmetries. First of all, from the fact that these 
equations do not have explicit x-dependence, they have the natural translation symmetry with 
respect to x — > x — a, where a is a real parameter. This means that if {P{t, x),Q{t, x)) is a 
solution for the equations, then {P{t, x + a), Q{t, x + a)) is also a solution. One more symmetry 
that naturally comes to our attention is the reflection symmetry with respect to re — > —x, which 
is a result of the invariance of the equations under the transformation d/dx — > —d/dx. This 
symmetry means that if {P{t, x), Q{t, x)) is a solution for the equations, then {P{t, —x),Q{t, —x)) 
is also a solution. Another feature that results from these symmetries is that (in the Hamiltonian 
picture) an initial data that is symmetric with respect tox^x + aorx^ —x maintains the 
same symmetry under the time evolution. This point will be discussed in the second subsection 
bellow. 

For the usual Gowdy model, both x ^ x — a and x —x transformations acting on 
{P{t,x),Q{t,x)) preserve the boundary conditions on P and Q, so that they express true sym- 
metries. For Nil and Sol models, however, these transformations do not preserve the boundary 
conditions, so that they do not express true symmetries. Nevertheless, we can find translation 
symmetries for Nil and Sol models in a generalized sense, and also find (generalized) reflection 
symmetries for them. We also find that the model admits larger amounts of reflection sym- 
metries than mentioned above. In the first subsection bellow, we obtain these translation and 
reflection symmetries in a systematic way. Based on this result, in the second subsection we 
discuss the dynamical point of view of the symmetries, which reveals a manifestation of the 
influence of spatial topologies to dynamics. 

4.1 Translation and Reflection Transformations 



All symmetries concerning the Gowdy equations ( p2D should have their origin in the metric (85). 
For example, this metric is invariant under the diffeomorphism x — > x — a, up to the redefinition 
of the metric functions {P{x), Q{x)) — >■ {P{x + a),Q{x + a)) (with the similar one for A). (For 
simplicity, we omit writing the argument t in the metric functions.) This redefinition accounts 
for the translation symmetry mentioned above. Recall that, if a metric is a solution for the 
vacuum Einstein equation, then so is the metric induced from a (spatial) diffeomorphism. Thus, 
in general, if such a diffeomorphism leaves the metric ( |85|) invariant up to a redefinition of the 
metric functions, this redefinition defines a symmetry for the Gowdy equations. 

Consider a diffeomorphism : {x,y,z) — > {(j)x{x,y, z), (l)y{x,y, z),(j)zix,y, z)) with t fixed. As 
remarked, (p niust preserve the characteristic form of the metric (^) so that a redefinition of the 
metric functions can make the metric invariant. In particular, to preserve the isothermal form 
e'^/^(— dt^ + dx^) of the metric we must have (f>x{x, y, z) = ±x — a. Moreover, {y, z) {(j)y, (j)z) 
must be a linear transformation on each y-z plane. That is, we must have 

(/): (x,y) ^ (±x-a,Ly), (42) 

where L £ SL^^) is a constant matrix and y is the column form vector of {y,z). Here, SL^'^^ = 
{m G GL(2,R)| det(m) = ±1}. (We do not have to consider the translations in the y-z plane 
of the type {x,y,z) (x,y + yQ,z + zq), which are isometrics for the metric (^), since these 
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give rise to no effect to the Gowdy equations.) We decompose the (jjx part into x ^ x — a and 
X —X, and will call symmetries involved with the former (generalized) translation symmetries 
and call symmetries involved with the latter (generalized) reflection symmetries. We can also 
consider symmetries generated by diffeomorphisms involved with the identity x ^ x, but since 
this type of symmetry is of little interest from the dynamical point of view discussed in the next 
subsection, we do not consider them here. We can find all possible translation and reflection 
symmetries defined above by finding (py^y^z) and 4)z{y,z) (or L) for each follows. 

We first consider the translation symmetries. As we will see in the next subsection, this 
symmetry implies that (in the Hamiltonian picture) a symmetric initial data with respect to 
a (generalized) translation maintains its symmetry under the time evolution. At the (locally) 
homogeneous limit, this implies that the diffeomorphism that gives rise to this symmetry should 
coincide with an isometry for the metric. Recall that there is a three-dimensional isometry group 
G for each case of T^, E'^, Nil, and Sol at the homogeneous limit. It contains a two-dimensional 
commutative isometry subgroup H ~ that is preserved as isometry group when the metric 
becomes £i^^-symmetric but inhomogeneous. Then, there is a one-parameter isometry subgroup 
K which acts transversely to H. In the £Z//^-symmetric but inhomogeneous cases, the actions 
of K are served as the desired translations. As easily found from the multiplication rules for the 
Bianchi I, VIIq, 11, and VIq groups, they are, respectively, given by 



Ta : {x,y,z) ^ {x - a, y, z), (43) 

Ea : (x, y, z) ^ (x — a, cos qa y -|- sinqa z, — sinqay -|- cos qa z), (44) 

Na : {x,y,z) ^ {x - a, y, z - ay), (45) 

Sa : {x,y,z) ^{x-a,e^i^y,e-^'^^z). (46) 

Since these leave the invariant one- forms for, respectively, Bianchi I, VIIq, II, and VIq types 
invariant, the corresponding transformations (redefinitions) takes the simplest form for the un- 
barred metric functions: 

Ta* ■ {P{x),Q{x)) (P(x -|- a),Q{x + a)) for metric functions in (85), (47) 

Ea* : {P{x),Q{x)) {P{x + a),Q{x + a)) for metric functions in (|3|), (48) 

Na* ■■ iP{x), Qix)) {P{x + a),Q{x + a)) for metric functions in (|l|), (49) 

Sa* : (P{x),Q{x)) (P{x + a),Q{x + a)) for metric functions in (|2^). (50) 

Here, we assume a = 27r/n for a positive integer n. Using the relations (|33|) , ([l9|) and (H), Ea*, 
Na* and Sa* expressed with the canonical metric (|85|) are found to be 

Ea* : {P{x),Q{x)) ^ {log^{P{x + a),Q{x + a),cosqa,-sinqa), 

^'(P(x -|- a),Q{x + a), cos qo, — sinqa) ), (51) 

Na* : {P{x),Q{x)) ^ {P{x + a),Q{x + a) + a), (52) 

Sa* : (P(x),Q(x)) ^ (P(x + a) -2qa,e2'i"Q(x + a)). (53) 



The functions <l> and ^ are defined in Eq. (p^) . Ta*, Ea*, Na* and Sa* preserve the boundary 
conditions for the £W^-symmetric models of, respectively, , E^, Nil and Sol types, so they 
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define the translation symmetry for each type. Note that the translations for each type form 
the infinite cyclic group Z for a fixed a, or form the cyclic group Z„ for a fixed a = 27r/n if the 
(periodic) boundary conditions are taken into account. 

Next, we consider the refiection symmetries, by which we mean the symmetries generated 
by the refiection transformations defined as follows. 

Definition Consider the map R : ix,y) — > {—x,Ly) for a given L £ SL^^\ The symmetry 
transformation R^: induced from R, acting on the variables (P, Q) of a given model, is called a 
reflection transformation, if the following three conditions are fulfilled; 

(i) For R (and i?*) to form the Z2-group together with the identity id, 

R^ = RoR = id. (54) 

(ii) As explained later, 

Conj{R)T = RoTo R-^ = T-\ (55) 

where T is the map inducing the translation transformation of the given model. (Hence, 
T is one of Eqs.(||) ~ (||).) 

(iii) i?* preserves the boundary conditions. That is, if variables (P, Q) satisfy the boundary 
conditions appropriate for the model, ii*(-P, Q) also satisfy the same boundary conditions. 

First, from condition (i), L must be one of 

_/l0\ ,_/-10\ ^ _ ( cos2d sin20 \ 

= 1^ 1 j ' = 1^ -1 j ' °^ = 1^ sin20 -COS20 ) ' ^^^^ 

where 6 G [0, vr) is a real parameter. These act on the y-z plane as, respectively, the identity, 
the rotation by angle vr, and the refiection with respect to the line Ig which passes through the 
origin with angle 9. Thus we have obtained the following candidates for the refiection R: 



Ro : {x,y,z) ^ {-x,y,z), (57) 
R77 ■■ {x,y,z) ^ {-x,-y,-z), (58) 
Rig : {x,y,z) ^ {—x,ycos26 + zsm26,ysm26 — zcos26). (59) 

Condition (ii) concerns a "compatibility" with the natural translation obtained earlier. That 
is, the refiection about a generic point x = a should be induced from that about x = by a 



translation. We can equivalently require that the conjugation, Conj(i?)T, of the translation T 
by a refiection R reverse the original translation, as indicated by Eq.(^). For example, the 
conjugation of the translation Ta for the model by the map Rq is given by Conj(i2o)?a = 
Rq o Ta o [Rq)~^ = Rq o Ta o Rq = T-a- This shows Rq is compatible with Ta- On the other 
hand, say, the conjugation of Na by Rq gives Conj(i?o)-^a = Ro ° Na o Rq = T_2a ° Na, which is 
not the reverse of A''^. All the compatible relations are given as follows: 

For r3 model: Conj(i?o )Ta = T_a, (60) 

Conj(P^)r, = r_„, (61) 
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For E"^ model: 
For Nil model: 



Coni{Rio)Ea = E^a (^G[0,7r))^ 
Coni{Rio)Na = N^a 



(62) 
(63) 
(64) 



For Sol model: Coni(i2,JS„ = 5_„ (e = -,-7r). (65) 

The transformations for the metric functions induced from these maps are given as follows: 

(66) 



Ro„R^* : iP{x),Q{x)) ^ iPi-x),Q{-x)), 

Ri^^ : (P(x),g(x)) ^ (log$(P(-x),Q(-x),-sin20,cos20), 

-^{P{-x), Q{-x), - sin 29, cos 29) ), 



where the functions <I> and ^ are defined in Eq. 
Sol {9 = 7r/4, 37r/4) models are 



The explicit forms for Nil 



Rio*, Rl^ 



/2* 



Rl^/4*, Rl-s^/4* 



(P(x),Q(x))^(P(-x),-Q(-x)), 



(67) 

0,7r/2) and 
(68) 



(P(x),0(x)) 



P + log(0^ + e- 



-2P\ 



Q 



Note that Rn 



to* = R-TT*, Riq* = RItt/2*^ and Ri^^^^ — ^•■'37r/4*' 
As a final task we are left with checking condition (iii). For the T"^ model, we at once see 
that if P{x) and Q{x) are periodic functions, then their images by i?o* and Ri^^ are also periodic 
functions. Thus, i?o* and i?;^* give reflections for the model. For the i?^. Nil and Sol models 
it is convenient to work with the unbarred variables. Interestingly the transformations (^) (for 
E^), ( |68| ) (for Nil) and ( ]69| ) (for Sol) are written with the unbarred variables in exactly the same 
form (except for the disappearance of bars); 
For E^: 



Rl, 



R 



la* 



For Nil: 
For Sol: 



R^ 



{P{x),Q{x)) {\og<^{P{-x),Q{-x),-sm29,cos29), 
-^{P{-x),Q{-x), - sin 26, cos 26) ), 



Rl,, : iP{x),Q{x)) ^ (p(-x),-Q(-x)), 

Q 



n/4* 



(P{x),Qix))^lP + logiQ' + e-'n, 



(70) 
(71) 
(72) 



From these it is trivial to see that the periodicity of the unbarred variables are preserved for E^ , 
Nil and Sol. This confirms that these transformations do give reflections for the corresponding 
models. Table || bellow summarizes the translation and reflection transformations for each 
model. 

Remark that we have inclusion relations. The model has the largest set of reflections. 
The E^ models have part of that for the model. The Nil and Sol models have parts of that 
for the E^ models. 
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Model 


Translation 


Reflections 




T 

a* 


Rq*, Rig* 


Ei 




Rig* 


Nil 


Na* 


Rio* 


Sol 







Table 1: The translation and reflection symmetry transformations. For reflections of and E^, 6 £ 
[0,^). 



4.2 Dynamical Interpretation 

As well known [Q], the Gowdy equations (^) admit a Hamiltonian formulation. Let 7rp{x) and 
t^q{x) be, respectively, conjugate momenta of P{x) and Q{x). The phase space V is spanned 
by the four functions {P{x),Q{x),7rp{x),TTQ{x)) (or similarly for variables with bars). Consider 
the set O of maps on "P, O = {o\o : V ^ V}. The Hamiltonian flow ipr, which generates the 
time evolution by time r, is naturally regarded as an element in O, V'r G O, for a flxed r. Vo is 
supposed to be the identity. 

What we are interested in here are maps (operators) which commute with t/j-r, so we define 

Ocom = {/ G OlV G R, foiP^ = ^p^of}. (73) 

We call an / G Ocom a symmetry operator. (If the operator / was a smooth flow like ipr we 
would be able to reformulate this commutativity to the vanishing of a Poisson bracket. However, 
since we consider discrete operators as / this is not the case. As a result, we do not actually 
need symplectic structures if only the flow tpr is properly defined.) For an initial data 7 S "P, 
iprij) is a solution for the Einstein equation as a function of r. If there is a symmetry operator 
/, then the function of r, /(ipril)), is also a solution, since this equals to tprifij)), which is 
the solution with the initial data f{j)- Therefore a symmetry operator generates a (in general, 
distinct) solution from a solution. 

Now, note that, while a symmetry transformation obtained in the previous subsection is 
defined in the configuration space spanned by P and Q, it induces an operator on the phase 
space V by differentiating the configuration variables with respect to r. We use in particular 
the Hamiltonian equations 

P = np, Q = e-^^TTQ. (74) 

(This relation is the same for all sets of barred or unbarred variables.) Thus, the symmetry 
transformations are also regarded as an operator on V. This operator is clearly a symmetry 
operator. Conversely, if there is a symmetry operator, it gives a symmetry transformation, as 
seen from the previous paragraph. Thus, in effect, symmetry transformations and operators are 
equal entities. Note that, since we have exhausted all the symmetry transformations concerning 
translation and reflection, we now know all the symmetry operators for them. 

To represent the symmetry operator obtained from a symmetry transformation we use the 
same symbol as that of the symmetry transformation. For example, the translation transforma- 
tion Ta* for model defined in Eq.(p7|) naturally gives the operator, denoted also as Ta*, that 
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simply shifts the argument by a constant: 

Ta* : {P{x),Q{x),Ttp{x),TtQ{x)) {P{x + a),Q{x + a),Ttp{x + a),TrQ{x + a)). (75) 

The translation symmetry operators Ea* (for E^), Na* (for Nil) and Sa* (for Sol) are also defined 
by similar maps if the phase space variables are taken as the "unbarred" ones. (cf. Eqs.(|4^), 
(@D, and dig).) 

An interesting consequence of the symmetry operators comes from considering data which 
are invariant with respect to this operator. For / G Ocom, let 

5(/)^{7GP|/(7)=7}. (76) 

For example, if / is a translation symmetry operator T2^/„*, the set 5(T2^/„*) is comprised 
of data in V which are translation-symmetric with period 2ii/n. If / is a reflection symmetry 
operator i?*, the set S{Rif) is comprised of reflection-symmetric data in V. The remarkable 
fact is that these symmetries are preserved dynamically, since Eq. (|73|) , together with 7(7) = 7, 
implies fi'ipri^f)) = i^ri^) for all r. In other words, any data 7 G S{f) remains within S{f) 
under the dynamical flow, i.e., '0r('5(/)) C S{f). Thus, S{f) defines an invariant subset of V 
for the flow ipr- 

We interpret the invariant subset S{f) as a dynamical character of the model. Motivated by 
this interpretation, we compare S{f) to see if there are distinctions in dynamical characters of 
the models. To this, recall that if we think of the phase space V as being spanned by the unbarred 
variables, V for every model is a space of four functions which are all periodic functions. (We 
think that the barred and unbarred variables coincide for the model.) Therefore the phase 
spaces for the four models are naturally identified. In this view, the flow ■0,- depends on the type 
of the model. Strictly speaking, however, since the "constraints" equations corresponding to 
Eq.([95|) take different forms for the E^, Nil and Sol models, this identification may be justified 
only approximately, but we neglect this point. 

First, consider the translations. As already mentioned the translation operator for every 
model is obtained by simply shifting the spatial coordinate like Eq.(^) for the unbarred vari- 
ables. Hence, all the invariant subsets S{Ta*), S{Ea*), S{Na*), and S{Sa*) for a = In/n, where 
n is a positive integer, are spanned by periodic data with period 2tt /n. In this sense, we roughly 
represent 

S{Ta,) ~ S{Ea.) ~ S{Na,) ~ S{Sa.). (77) 

Therefore we conclude that for the property of preservation of translation symmetry, there is no 
significant distinction between the four models. 

Our main interest is therefore in the reflections. As we remarked, reflection operator 

i?o* : {P{x), Q(x), 7fp(x), 7fQ(x)) ^ {P{-x), Q(-x), 7fp(-a;), 7fQ(-x)) (78) 

does not preserve the boundary conditions for E^ , Nil and Sol models, so that it is not relevant 
for these three models. At first sight, the similar operator acting on unbarred variables 

i?^, : {P{x),Q{x),T:p{x),nQ{x))^{P{-x),Q{-x),np{-x),T:Q{-x)), (79) 
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which does preserve the boundary conditions for the three models, seems to be a natural reflec- 
tion operator for them. Indeed, we can regard the spatial coordinate x for the unbarred variables 
as a natural coordinate in view of the translation property shown above. That is, the reflection 
with respect to this coordinate seems like the most natural one. Nevertheless, the reflection 
symmetry with respect to this operator imposed on an initial data for E^, Nil or Sol model 
is not preserved under the time evolution. This fact can be seen directly from the dynamical 
equations for the unbarred variables. For example, those for the Nil model are given by (See 
Eq-dl)) 

Q - e-^^Q" + 2{PQ - e'^^P'iQ' - 1)) = 0, (80) 

which are apparently, in contrast to the model, not invariant under the transformation 
d/dx — > —d/dx, due to the factor {Q' — 1). Similarly, those for the and Sol models are not 
invariant under the same transformation. As a result, in contrast to translation, there are no 
invariant subsets for reflection in the E^, Nil and Sol models which naturally correspond to the 
invariant subset 5(i?o*) for the T'^ model. 

However, the model admits another one-parameter family of nontrivial reflection oper- 
ators Rig*. (See Eq.(|67|)) The momentum part is derived from Eq.(|67|) by differentiating with 
respect to r and using Eqs.(^4[). The result is given by 



Rie* ■ (7rp(x),7rQ(x)) 



((c + QsY - e-^^52)7fp + 2e-'^s{c + Qs)7rQ 
(c + Qs)'^ + e-'^Ps^ 

2s{c + Qs)Ttp - ((c + Qs f - e~^^s'^)7tQ 



(81) 



where c = cos 29, s = sin 29. 

As shown in the previous subsection, for Nil and Sol models only isolated values (0 = or 
7r/2 for Nil, and 7r/4 or 37r/4 for Sol) of 9 are permissible. As a result, the invariant subsets for 
reflection in the Nil and Sol models are much limited compared to the model. More precisely, 
we have the following. 

Proposition 9 Let S{Rt), S{Re), S{Rm), and S{Rs) he the unions of the all invariant subsets 
for reflection in, respectively, the , E^ , Nil, and Sol models. Then, the inclusion relation 

S{Rt) D S{Re) D iS{RN), SiRs)) (82) 



holds (if neglecting the constraints corresponding to Eq.(95)). 



Proof This can be seen from S{Rt) = 5(i?o*) U {UgS{RigJ), S{Re) = UeS{RigJ, S{Rn) = 
S{RiJ, and S{Re)=S{Ri^^ J. □ 

It should be stressed that the relation ( |8^ ) truly reflects the dynamical properties of the 
models. Note that, since the (periodic) boundary conditions imposed on initial data are the 
same for all models if the phase space V is spanned by the unbarred variables, a map / which 
preserves the boundary conditions for a model always preserves the boundary conditions for 
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other models, too. This means that there are the same amounts of "reflection symmetric initial 
data" for the four models. However, the time evolution does not always preserve the symmetry 
of such an initial data, as we have seen. To conclude, we have clarified how much "reflection" 
symmetric data exist for which the symmetry is preserved under the time evolution for each 
model, and found, in particular, the inclusion relation (|8^). This may be interpreted as a 
manifestation of the infiuence of topology to dynamics. 



5 Summary and Comments 

We have made the structures of £Z^^-manifolds clear and classified the possible topologies of 
them. For convenience we have split the set of these manifolds into two kinds; the first kind, 
those with local Killing vectors which are not degenerate everywhere, and the second kind, those 
with ones which are degenerate somewhere. The local Killing vectors of any £W^-manifold of the 
second kind are actually defined globally, so that all the £Z//^-symmetric models of the second 
kind are contained in the usual Gowdy models. On the other hand, it is only that is contained 
in the Gowdy models among the varieties of the /Z/^-manifolds of the first kind, so that in this 
paper we have basically restricted ourselves to the models of the first kind. 

In the case that three-manifold M is an £Z//'^-symmetric space of the first kind, M is a 
T^-bundle over the (Theorem ^), and according to the corresponding geometric structure, 
M is naturally characterized by one of E^, Nil, or Sol. In each case, the possible topologies 
of M can be classified more precisely, but if we are interested in the dynamical properties of 
the corresponding spacetime models it is not necessary to consider all the models one by one. 
Restricted to representatives for this "dynamical equivalence class," the number of models to 
consider is to great extent decreased. In particular, there is three representatives for E"^ (cf. 
Sec. 3.3), and there is only one for Nil (cf. Sec. 3.1). The representatives for Sol are parameterized 



by only one discrete parameter q (cf. Sec. |3.2D . 

We have given two ways of representing the metric for each of E^, Nil and Sol models. (We 
distinguish between the E^ and models. See Convention in Sec. |3.3| .) Note that the metrics 
of all the £Z//^-symmetric spacetimes can be represented locally in the same canonical form, but 
with distinct boundary conditions. In this view, however, global symmetries that arise from 
the spatial topology tend to be unclear. On the other hand, the other way of representing the 
metric makes the geometric structure the i2Z//^-manifold M admits manifest. Geometrically it 
is obtained by "relaxing" the corresponding locally homogeneous metric. This type of metric 
has another advantage that it gives rise to a natural reduction of the model; that is, since the 
boundary conditions are always given by the periodicity for the spatial coordinate, the spatial 
manifold naturally reduces to the S^, the base space of the bundle. Note that this fact makes 
the identifications of the phase spaces for the varieties of the models possible. 

Finally, as an application of these fundamental facts we have given the translation and re- 
flection operators which commute with the time-evolution, and have discussed their significance. 
Since, as stated above, the metrics of £Z//^-symmetric models can always be locally represented 
in the same form, these models are considered to have the same local dynamical properties. 
However differences of topologies affect the global dynamical properties of the models through 
the differences of the boundary conditions. In this paper we have examined whether or not a 
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global symmetry imposed on an initial data is preserved in time, or how much there exist global 
symmetries which are preserved in time, for each model of the first kind. As a result, we have 
indeed found that there are remarkable distinctions in the properties of reflection. The freedom 
of reflection symmetries is the largest as for the model, and we have the inclusion relation 
(|8^). In particular, naive (even type) reflection symmetric initial data for the E^, Nil or Sol 
model do not evolve maintaining its symmetry in time, in contrast to the model. 

In the following we make some comments. A first one concerns another point concerning the 
correlations between topology and dynamics. As for the dynamics of the Gowdy model, what 
has been being taken an interest most is whether the AVTD conjecture jl^ is correct or not 
that predicts a universal behavior of the approach toward the initial singularity. This conjecture 



has been basically supported both analytically 2C] and numerically pi|, p^. However, the 



subtlety has been pointed out [21| that there is a measure-zero set of spatial points where the 
AVTD behavior is not achieved. Here we concern ourselves with these nongeneric points. It is 
known that at these points Q' = 0. Since this condition is local, it does not depend on the 
topology. However, note that while points of Q' = are inevitable in the model, since Q{x) 
is a (smooth) periodic function, such points are not necessary in the Nil and Sol models. Hence 
we can naturally expect that topology affects the tendency of the appearance of the nongeneric 



points. The points where Q' = correspond, from Eqs.(f[9D and (|29|), to the points where 

Nil: Q' = 1, (83) 

Sob (log|Q|)' = 2q, (84) 

for the unbarred variables. The unbarred Q{x) must be periodic, but this condition does not 
force the existence of points such that the above conditions are fulfilled. However, since the 
evolution of Q freezes when approaching the initial singularity, we can naturally expect that 
nongeneric points are still generated if the maximal gradient of Q at an initial surface is large 
enough and if, as a result, Q^ax > 1 (^il); (loglQI)max > (Sol) at a time when the 
evolution of Q freezes. Indeed, this property has been observed from numerical simulations the 
author performed. A detailed description for this point will be reported elsewhere. In Ref. |^ it 
was suggested that the nongeneric points are not generated in the Nil and Sol models, since the 
conditions corresponding to Eqs.(p^) and ( |8^ ) were written in such different forms with different 
choice of metric functions that the conditions seem not to be fulfilled. This claim is not correct, 
as seen from above. Still, we may claim that the -E^ model is most likely to generate nongeneric 
points. 



Next we comment on the locally U(l)-symmetric models |23| that are less symmetric than 
the /^Z^^-symmetric models. This model admits only one spatial local Killing vector. The spatial 
manifold M is again assumed to be closed. First, from an analogous analysis to the proof of 
Lemma the Killing orbits are found to be closed and are homeomorphic to S^, if the local 
Killing vector does not degenerate everywhere (that is, if the model is the "first kind"). Hence 



M is naturally a Seifert fiber space [17|, and still admits a geometric structure, which is one 
of E^, Nil, X R, SL(2,R), 5^ x R, or 0. In the case that M is the "second kind", 
i.e., the local Killing vector vanishes on somewhere, the manifold M is obtained from a Seifert 
fiber space by performing a finite number of particular kind of Dehn surgeries. The resulting 
manifold M is found (from one in the series of famous theorems of Thurston concerning the 
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structure) to again admit a geometric structure. Detailed accounts and applications will be 
presented elsewhere. 

A Vacuum Einstein equations for Gowdy Spacetimes 

In this Appendix, we summarize the standard prescription |^ of the reduced vacuum Einstein 
equations for Gowdy spacetimes, together with generalizations to Nil and Sol types. 

We consider the two-surface orthogonal metric (|l^). This is used most frequently in the 
literature since the Einstein equations become to great extent simpler. 

More explicitly, we write 

ds^ = e^/^{-dt^ + dx^) + R{e^{dy + Qdz)^ + e'^dz'^). (85) 

where P and Q are functions of t and x. Then, the vacuum Einstein equations for the metric 
functions are explicitly given as follows. The dynamical equations are 

Ptt - P" -e^^{Q^t-Q'^)+R-\RtPt-R'P') = 0, (86) 
Qtt - Q" + 2{PtQt-P'Q') + R'\RtQt-R'Q') = 0, 

and 

Rtt - R" = 0, (87) 

where subscript t represents the derivative with respect to t and the primes represent the deriva- 
tives with respect to x. The remaining constraint equations takes a simple form in the null 
coordinates u = t — x and v = t + x; 

R{Pl + e^^Ql) = -2Ruu + Ruilu + RT^Ru), 

R{P^ + e^^Ql) = -2R^^ + R^{% + R-^R^), (88) 

where the subscripts stand for the derivatives thereof. 

The last two equations can be solved for 7^ and % at any spacetime point wherever R^ 
and R^ do not vanish. In this case the integrability condition for 7, namely 'juv — Ivu = 0, 
is automatically satisfied if the dynamical equations (^) are satisfied. Hence we can integrate 
the equations and obtain 7 by using solutions of the (unconstrained) dynamical equations (p^). 
If there is a spacetime point such that ^„ = (or = 0), Eqs.(^) cannot be solved for 7^ 
(or 7t,), so that the constraint equation ( p^ ) constrains P and Q at that spacetime point. This 
condition is called the "matching condition" 0. Explicitly, 

Pl + e^^Ql = -2R-^Ruu at i2„ = 0, 

Pv + e^^Ql = -2R-^Rvv at Ry = 0. (89) 

From the observation that the left hand sides of Eqs.(p9|) are positive semidefinite, we find 
that for E > 0, 

Ruu < at Ru = 0, 

< at R^ = 0. (90) 
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An implication of these inequalities is called the "corner theorem" Here, a corner is a point 
in the t-x plane such that Ru = ov = 0, that is, a point where the gradient of a level curve 
of R becomes null. The implication of the inequalities would be clear if observing some level 
curves of R near a corner in the u-v coordinates. It should be stressed that these inequalities do 
not depend on the boundary conditions for P, Q, and 7, so hold for any spatial topology. (This 
theorem, however, tacitly assume the finiteness of 7^ and % at the corner. If we allow for Ruju 
and Rvjv to remain finite at the corner with diverging 7^ and 'jy, we may have counterexamples 
to this theorem. p5[) 

For any i2Z//^-symmetric spacetime of the first kind, the area function R must be periodic 
(cf. Eq.([l3D). Hence, corners appear only in pair, implying that R cannot have a corner. We 
can therefore choose R spatially constant 

R{t,x)=t, (91) 
(using the remaining freedom of the coordinate transformations of type u F{u),v G{v)). 



In this case, the reduced Einstein equations (86) and (|8q ) become 



with 



p _ - e2^(0' - e-2-0") = 0, (92) 

Q - + 2(P0 - e-2"P'Q') = 0, 



p\e-^^P" + e'^^{Q +e-^^Q'^), (93) 

2P, 



X' = 2{P'P + e^^Q'Q). 



Here, we have put 



^ = e-^ ^ (94) 

Dots in these equations represent derivatives with respect to r (not t). A should be periodic for 
any £Z//^-symmetric spacetime of the first kind. This leads to the constraint for P and Q, 

{P'P + e^^Q'Q)dx = 0, (95) 





because of the second equation of Eqs.(p^). 

The functions P and Q are not periodic functions for the Nil and Sol models. Appropriate 
boundary conditions are obtained from Eqs.(|l9|) and (|29|). Specifically, 
For Nil: 

P{x + 27r) = P{x), Q{x + 27r) = Q{x) - 2-k. (96) 

For Sol: 

P{x + 2^) = P{x) + 4^q, Q{x + 27r) = e-^^'iQ(x). (97) 

Remark The last boundary conditions (p^) and ( |97| ) do not necessarily coincide with 
those obtained from Eqs.([T^). This is because the choice of the covering group T (cf. Eq.([lO|)) 
is different. In obtaining Eqs.(p!^) we have fixed T first and found the appropriate metric, but 
for Eqs.(^) and (|97| ) we try to fix the universal covering metric with appropriate F as long as 
possible. This is the essential point in the dynamical equivalences for the Nil and Sol models 
presented in Secjs]. 
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B Compactification of Sol 



In this Appendix we show calculations needed to perform the compactification of Sol. This is 
basically a review of that presented in Ref. |Q . In this reference, the form of mcg_,_ (T^) ~ SL(2, Z) 
was restricted to that of Eq.@, and the parameter q (See bellow) was not introduced. In this 
review we show an explicit calculation with general form of the matrix and with q, though no 
essential difference appears. 

What we do is to embed the fundamental groups into Sol~ Gvi,, (the Bianchi VIq group). 
For notational convenience, letting a = g2, b = gs, and c = 51 in the representation (^), 



TTl 



a,b,c; [a,b] = l,cac =a^b^,cbc = a'^b'^ 



where 



A = 



p q 

r s 



eSL(2,Z), with |Tr^| > 2. 



(98) 
(99) 



We represent the tti -generators with their components in Gyio 



ai \ 








, b = 




03 / 




[ 63 




(100) 



To avoid confusions with powers, we use subscripts to distinguish each components as above. 
The multiplication rule for GvIq is given by 



a \ 




X 


?) 




y 




v 


z 



a + X 
7 + e<5"z 





( —a 











(101) 



where q > is a parameter (, which has nothing to do with the elements in the matrix A). Left 
invariant one- forms are then given by Eq.(p6|). 

First, note the first components of the relations 



cac 



aW, cbc-^ = a%'. 



(102) 



Noting that the first component of the product of two elements of Sol is simply given by the 
sum of the first components of the elements, we obtain oi = pai + rbi and 61 = qai + sbi, that 
is. 



p-1 



r 

s - 1 



ai 
bi 



(103) 



The determinant for the matrix appearing in the left hand side is given by det A — TrA = 1 — Trvl, 
which is less than —1 (when TrA > 2) or greater than 3 (when TrA < —2), from the assumptions 
on the matrix A. Hence, the matrix has its inverse, implying 



ai 



0. 



(104) 
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Under this condition the relation [a, b] = 1 becomes trivial. Moreover, we can at once calculate 





coc 





a2e~'^^^ I , etc., 



so that the second and third components of the relations ( 102 ) are given by 



(105) 



A 



02 



pqci 



02 

b2 



A 



as 



03 



(106) 



These reveal that ( | and ( ^'^ | are eigenvectors of A, and e^*^*^^ are the eigenvalues. The 

\b2 J \b3 J 

characteristic polynomial is — Tr^ A + det A = 0, or equivalently. 



A^-TrAA + 1 = 0. 



(107) 



Since the eigenvalues e^^^^ must be positive, the embedding is found to be possible only for 
the case Tr^ > 2. However, even for the case TiA < —2, if the map 



h : ix,y,z) {x,-y,-z) 



(108) 



ci 

can be regarded as an isometry, we can perform the embedding by putting c = /lo | C2 | , where 

the column vector in the right hand side is an element of Sol, and o represents the composition 
of maps. It is an easy task to retrace the calculation above. We will find in particular that the 
eigenvalues of A must equal to —e^^'^^ rather than e^^'^^ . 
Finally, the formal solutions for embedding follows. 

(i) Case of TrA > 2: Letting the eigensystems (the pairs of an eigenvalue and the corre- 
sponding normalized eigenvector) for the matrix A be 



pqci 



U2 

Us 



-qci 



V2 
VS 



(109) 



the solution is 






(110) 



where /?, 7, C2, and C3 are parameters. Note that, while the product qci is directly connected 
with the eigenvalues of A, we can freely specify one of g or ci. 
(ii) Case of TiA < —2: Letting the eigensystems be 



U2 

Us 



~qci 



V2 

vs 



(111) 
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( PU2 \ 


( 


a = 




, b = \ -fV3 




\ ) 


V 




the solution is 

f [3U2 \ ( (3U3 ^ 

/lo I C2 I , (112) 

where /3, 7, C2 and C3 are parameters. As in the case (i), one of q or ci can be freely specified. 

Remark We have seen that the fundamental group for Tr^ > 2 can be embedded into 
Gvio- Also, we have seen that the fundamental group for TrA < —2 can be embedded into G-yj^, 
the group generated by Gyio and the Z2-factor h. It is easy to see that the same groups can 

be embedded into the smaller group -ffyio C Gyio which is given by making the first component 

(2) (2) (2) 

of Gvio discrete like ciZ, or into -f^vio' where -f^vio ^vio group generated by Hyig and 



(2) 

the map h. Since -ffyio isometry group of the universal covering spacetime (|27D (with 



the appropriate choice of q, as shown in Sec.^), the same embeddings are possible for this 
(inhomogeneous) spacetime. 
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